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Abstract: We study the phenomenology of a light stop NLSP in the presence of large mixing 



j^ with either the first or the second generation. R-symmetric models provide a prime setting 

for this scenario, but our discussion also applies to the MSSM when a significant amount of 
mixing can be accommodated. In our framework the dominant stop decay is through the 
flavor violating mode into a light jet and the LSP in an extended region of parameter space. 
There are currently no limits from ATLAS and CMS in this region. We emulate shape-based 
hadronic SUSY searches for this topology, and find that they have potential sensitivity. If the 
£<■} extension of these analyses to this region is robust, we find that these searches can set strong 

■^ exclusion limits on light stops. If not, then the flavor violating decay mode is challenging 

and may represent a blind spot in stop searches even at 13 TeV. Thus, an experimental 
investigation of this scenario is well motivated. 
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1 Introduction 

With the discovery of the particle consistent with standard model (SM) Higgs boson [1, 2], the 
hierarchy problem takes on a concrete and immediate nature. Low energy supersymmetry is 
one of the most studied and robust solutions to the hierarchy problem. A simple realization of 
supersymmetry is the Minimal Supersymmetric Standard Model (MSSM), which in its most 
general form suffers from a very acute flavor problem [3, 4]. The tension between flavor and 
natural low energy SUSY can be economically solved in models which predict spectra with 
degenerate squarks (and sleptons). In the pre-LHC era this paradigm provided a compelling 
framework of flavor in low energy SUSY models. However, the LHC bounds on these models 
push the spectrum of new states to the TeV scale, already requiring the SM to suffer some 
amount of fine-tuning. 

Natural SUSY models surviving the LHC bounds now require a non-trivial flavor struc- 
ture for the squark mass matrix so as to produce a hierarchical spectrum with the first and 
the second generation significantly heavier than the third generation [5, 6]. In the MSSM the 
structure of possible squark mass matrices is severely restricted by low energy observables 
[7-10]. R-symmetric models provide a compelling setting for investigating more complicated 
flavor structures since they generically admit larger mixing among the squarks than in the 
MSSM [11]. This is a consequence of the fact that most flavor violating observables turn 
out to be proportional to terms which break R-symmetry, and are hence suppressed [11-13]. 



Therefore, a significant amount of flavor mixing is possible and quite generic. This could 
be an interesting ingredient in building flavorful SUSY breaking/mediation models, and also 
significantly alter the LHC signatures of SUSY particles [14]. 

Models with an (approximate) R-symmetry [11, 15-18] are additionally interesting as 
an alternative to the MSSM since they require the gauginos to be Dirac fermions.* This 
makes them less minimal than the MSSM in terms of the particle content, adding adjoint 
superfields for each SM gauge group, but on the other hand the Dirac character of the 
gauginos (particularly gluinos) can significantly soften the LHC exclusion limits on direct 
squark production [19-21]. Moreover, Dirac gluinos above the TeV scale are natural [22], 
while Majorana gluinos above that scale begin to saturate the naturalness threshold [5, 23]. 
This is a particularly interesting feature in light of the LHC limits on gluinos, ~ 1.5 TeV [24]. 

In this paper we will investigate the implications of large flavor mixing for the LHC 
stop phenomenology. R-symmetric models provide an optimal setting for this scenario, 
but our discussion also applies to the MSSM when a significant amount of mixing can be 
accommodated. We will focus on a light stop NLSP (Next-to-Lightest Supersymmetric 
Particle) with a mass splitting AM = mj — misp ^ 250 GeV with the LSP. In this region 
of parameter space the mixing induced decays t —> j+LSP dominate, and hence completely 
alter the phenomenology of the stop. 

The paper is organized as follows: we briefly review the main features of models with an 
approximate R-symmetry. We then present the effect of large flavor mixing on squark decays, 
focusing on the stop. The rest of the paper is devoted to the collider phenomenology of this 
scenario, with emphasis on the region of parameter space where the flavor violating decays 
dominate. 

2 Mixed stops in R-symmetric models 

One of the main features of R-symmetric models is that gauginos are Dirac fermions instead 
of Majorana. The R-symmetry also forbids A terms, which implies that there is no left- 
right mixing. Further, the //-term is also not allowed. The combination of these features 
relaxes the SUSY flavor problem, allowing for a larger flavor violation both in the squark 
and in slepton sector than in the MSSM (see [11] and [13]). For instance, models with Dirac 
gluinos contribute to meson mixing through dimension 6 operators which are hence suppressed 
relative to the dimension 5 operators present in models with Majorana gluinos allowing for 
a potentially large mixing. The imaginary part of the K — K mixing is severely constrained 
by the measurement of ex- Therefore, the mixing in the first and second generation requires 
mild flavor structure in order to be consistent with the limits. 

Flavor changing processes like fi — > e^f or b — > 57 require a helicity flip in the correspond- 
ing diagrams which typically would arise from either a Majorana mass insertion, left right 
mixing between the squarks or a /i-term insertion in the large tan/3 limit; all these terms are 

* The R-symmetry cannot be an exact symmetry since it is at least broken by the gravitino mass. The 
effects of the R-breaking are therefore communicated to the visible sector through anomaly mediation, see [11]. 



forbidden by the R-symmetry. Therefore, one is left with a mass insertion in the external line 
which is significantly smaller, leading to the suppression of all the AF = 1 processes. The R- 
symmetry also protects from dangerous contributions to AF = 1 processes like b — > 57 arising 
from a light stop. Moreover, the dangerous one- loop contributions to EDMs are forbidden by 
the R-symmetry since they all require either a Majorana mass insertion or squark left-right 
mixing. 

The R-symmetry does not completely solve the flavor puzzle in SUSY. A mild flavor 
structure (in particular, a small phase) is still needed in order to satisfy the bounds from 
K — K mixing. However, many flavor violating observables are absent, and others signifi- 
cantly suppressed. Importantly, these suppressions are obtained independent of the SUSY 
breaking/mediation paradigm (beyond the requirement of R-symmetry) . Therefore, the set of 
flavor assumptions for a viable R-symmetric SUSY model is much smaller than in the MSSM. 

In this work we are interested in exploring the consequences of having large flavor mixing 
on the stop phenomenology. As a simplifying assumption we only consider mixing between 
the second or the first generation and the third generation squarks, both in the left handed 
and in the right handed sector. We assume the unmixed generation to be heavy. Qualitatively 
our results continue to hold for more general flavor structures. 

The state responsible for canceling the top quadratic divergence is the "33" state in 
the squark mass matrix M~. Therefore, naturalness requirements affect (M~)s3, which we 
choose to be < (400GeV) 2 . In presence of flavor mixing (M?) 33 is not the physical mass of 
the (mostly) stop squark, but is related to the two mixed squarks masses by the following 
relation: 

(M?) 33 = cos 6 2 m 2 + sin 6 2 m 2 (2.1) 

where j is either a (mostly) charm or up squark. The Higgs mass parameter thus becomes 
sensitive to the potentially heavier squark mass scale through the mixing. Therefore, the stop 
mass is not a robust measure of fine tuning in this scenario, and the fine tuning is in general 
worse than that deduced from the stop mass. The estimation of fine tuning then involves 
both the mixed generations. This is beyond the scope of this paper, and in the following 
we will assume the other generations to be above 650 GeV making them safe from the LHC 
bounds on three degenerate squarks (flavor-mixed partners of tn and ti). This limit can be 
derived by rescaling the bounds on 8 degenerate squarks with decoupled gluinos. As shown 
in figure 1, this spectrum is still compatible with naturalness {(M"-)^ < (400 GeV) 2 ). 

It is also important to note that the presence of a large mixing angle can potentially 
affect the production cross section of stops at the LHC, since they can be now produced via 
their mixing with the charm or up-quark. However, for sufficiently heavy (Dirac) gluinos we 
can assume this to be a subdominant effect, especially for the case when the stop mixes with 
the charm. Having heavy gluinos also means that the only source of squarks production at 
the LHC (at 7-8 TeV) is direct q*q production. 
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Figure 1: Natural region, which is chosen to be correspond to (Mq)^ <400 GeV, for m~- 
650 GeV (light purple) and m~- = 1 TeV (darker purple). 



3 Flavor violating stop decay modes 

In this section we discuss the impact of large squark mixing in the third generation phe- 
nomenology, focusing on the NLSP stop scenario. This first requires us to specify the nature 
of the LSP. Third generation final states are preferred if the LSP contains an appreciable 
Higgsino component. This means that the flavor violating decays dominate for the first and 
second generation squarks, while the third generation phenomenology is not affected by the 
mixing. In this work, we focus on the case where the coupling is flavor universal, and the 
LSP is a pure bino neutralino or a gravitino.^ The gravitino is often the LSP in R-symmetric 
models since its mass is related to the amount of R-breaking communicated to the visible 
sector through anomaly mediation [11]. Alternatively, a pure bino LSP is a good dark matter 
candidate. It is a pseudo-Dirac fermion, which means it annihilates more efficiently than the 
Majorana bino leading to the correct relic density in a large region of the parameter space 
[25]. On the other hand at the present day the pseudo-Dirac bino behaves like a Majorana 
fermion making a it safe from direct detection constraints [26, 27]. We now discuss the 
phenomenology of a light stop m^ < 400 GeV for these two scenarios. For this purpose it is 
useful to separately discuss two regions: m^ > m t + ra^sp and m 4 - < m t + rriLSp. 

*We do not consider a wino LSP since it typically implies light charginos, which are not the focus of our 
study. 



3.1 m t ~ > m t + m LS p 

In this region of the parameter space a stop decays either into a t+LSP or j+ LSP where 
the jet is either a charm or up type jet (for large #23 or #13 respectively, denoted 9 in the 
following). For a sufficiently heavy stop the branching ratio for the decay mode into light jet 
is sin 2 9, but in the almost degenerate region m^ ~ mt + mLSP it gets enhanced due to the 
phase space suppression of the decay into top quarks, and hence could become the dominant 
decay mode. For instance, for a massless bino LSP the branching ratio for these decay modes 
are: 

~ sin 9 
Br(i^jB) = — ^-, (3.1) 

sin 2 + cos 2 0(1 -^|) 2 

t 

cos 2 9(1 -=|) 2 

Br(t^tB) = — *— 5-, (3.2) 

sin 2 9 + cos 2 9(1 -^) 2 

i 

and for a massless gravitino LSP: 



sin 2 ■ 



Br(t^jG) = — ^-, (3.3) 

sin 2 9 + cos 2 9(1 - ^f) 4 

cos 2 #(l-4) 4 

Br (t^tG) = — m ~t . (3.4) 

sin 2 9 + cos 2 9(1 - ^|) 4 

t 

These equations apply for decays of both tj_, and £#. 

As shown in figure 2 the FV decay mode dominates over t+LSP for mass splitting AM ~ 
mt due to the large mixing angle and phase space suppression. For a mixing of sin 9 = 0.25 
the FV has order one branching ratio for stop masses up to 250 GeV for a massless bino LSP. 
For the gravitino this mode continues to dominate for much smaller mixing angles. Therefore, 
the FV mode is significant up to masses where dedicated searches for the topology £+LSP 
starts to be sensitive. 

This is interesting since the region mt < m^ < 250 GeV represents an important lack 
of coverage of LHC analyses. Indeed the sensitivity of all the stop dedicated searches in the 
channel t plus LSP starts at 230-250 GeV both for ATLAS [28] and CMS [29] for a massless 
LSP. The challenge here is in distinguishing a possible stop signal over the top background. 
Due to the experimental difficulties this region (mj ~ mt) is called the stealth stop region 
[30]. For proposed dedicated searches in this region see [31, 32]. Alternatively, one can hope 
to discover the stealth stop relying on a different decay mode. In the MSSM when the stop is 
the NLSP the only other possible decay mode is into an off shell top and the LSP. However, 
the branching ratio for this three-body decay ends abruptly at m ( - = mt for a neutralino LSP. 
If the gravitino is the LSP instead the branching ratio for three-body decay is significant 
(~ 20%) up to 200 GeV, so the stealthy region is less broad [30]. We have shown that in our 
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Figure 2: Contour plot for the branching ratio Br(t — > B + j) (solid) and Br(t — > G + j) 
(dashed) as function of the stop mass and the mixing. Contours are shown for branching 
ratio 75 % (blue) and 50% (red). The LSP is taken to be massless in both cases. 

scenario, if the mixing is significant, a stop almost degenerate with the top quark will decay 
mostly into a light jet plus LSP. This can simplify the challenge of looking for stops into the 
stealth region. We will discuss this in the next section. 

For heavier stops, m^ 3> mt + tulsp, the branching ratio for the decay mode into 
light jet is still large, but subdominant (it is proportional to sin 2 #). As discussed in [33] 
this can potentially relax the bounds on the stops from the t LSP searches, [28] and [29]. 
Another interesting phenomenological implication of a significant squark mixing is single top 
production at the LHC (pp — > t*t with t —> j+LSP and t* — > t*+ LSP and vice versa). This 
topology constitutes a "smoking gun" signature for scenarios with a large squark mixing and 
it has been discussed in [33] and more extensively in [14]. It is important to note that the 
CMS and ATLAS analysis do not have a dedicated analysis for this mixed decay topology 
despite the low SM background. 

3.2 m t ~ < m t + rriLSP 

The two-body flavor conserving stop decays are forbidden for m^ < mt + rriLSP- We will show 
that, for sufficiently large mixing angle, the two-body flavor violating mode dominates over 
the three-body decay i^-bW LSP. 
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Figure 3: Contour plot for 90% branching ratio for the two-body decay ti — > jLSP for a 
bino LSP (solid) and a gravitino LSP (dashed), for m^- — rriLSP = 150 GeV. A similar plot 
applies to the tji decay modes. 



If the gravitino is the LSP then the decay width for two open decay modes are [34] : 

a {mi — mw) 



T(t -»• bWG) ~ cos^ 6 



r(t^jG)~sin 2 



sin 2 6 W 1287r 2 m2 F 2 ' 



W J 



TO- 



16ttF 2: 



(3.5) 
(3.6) 



where VF is the SUSY breaking scale. The branching ratio for the flavor violating (FV) 
two-body decay is: 

1 



Br{t -* jG) 



i + cotg 2 ( s ! 2fl m s 



(3.7) 



When the gravitino is the LSP the two-body decay mode dominates over three-body 
decays even for very small mixing angles, as shown in figure 3. This decay mode is prompt 
for a low SUSY breaking scale {yF < 100 TeV). For these SUSY breaking scales the gravitino 
is extremely light and hence in the following we will take it to be massless. If the LSP is 
purely bino, the two-body decay dominates (with branching ratio greater than 90%) over the 
three-body decay for a larger mixing angle (around sin 2 9 ~ 0.08) than in the gravitino case 



as it shown in figure 3 for a left handed stop (the same conclusion applies to tp) . We conclude 
then that the FV two-body decay mode is the dominant decay mode for the stop in the entire 
region of the parameter space where mj < mj + m-isp for moderate mixing. Also, as shown 
above, the region where it dominates extends above the m^ = rrit + ttilsp threshold into the 
so-called stealth region. Therefore, the decay mode into jet plus LSP is a potential discovery 
mode for a stop NLSP almost degenerate with the LSP. It is interesting to briefly compare the 
FV stop decay mode in our scenario with the MSSM with minimal flavor violation (MFV). 
In the MSSM the dominant decay mode is the 3-body decay into bW plus the LSP, and when 
this is kinematically closed (m^ — ttilsp < rnw + m b) the dominant decay mode is either the 
loop suppressed two-body flavor violating (FV) decay mode, t — > c+LSP or the four body 
decay. The FV stop decay mode is suppressed compared to the three-body decay since the 
effective coupling t xl c is of order 10 -5 assuming MFV (see [35-37]). 

4 LHC sensitivity to a light stop NLSP with large flavor mixing 

In the previous section we showed that the FV decay mode into a light jet (either a charm or 
up-quark initiated) and the LSP is the dominant decay mode and hence a potential discovery 
mode for a stop almost degenerate with the LSP when a moderate squark flavor mixing is 
present. It is therefore interesting to determine the LHC sensitivity for this topology, and the 
bounds on the stop mass in this region of the parameter space. 

The Tevatron already puts interesting constraints on stops decaying into charm jets. 
Both DO [38] and CDF [39] have performed a dedicated analysis for t — > c+LSP, putting a 
bound on the stop mass as strong as 180 GeV for a 90 GeV LSP. However, both the more 
compressed region AM < 40 GeV and the region outside m^ — m^ < my/ + mi,, were left 
unconstrained. These analyses rely on c-tagging, and hence do not directly constrain the FV 
stop decay into an up quark. 

ATLAS and CMS do not have any dedicated analysis for this topology — light stops are 
searched in the \\ b channel [40]. Nonetheless, since our final state contains missing energy 
and two jets, jets plus missing energy (MET) LHC searches are potentially sensitive to this 
topology. Limits on our case can be derived from simplified models for squark searches with 
decoupled gluinos. However, CMS and ATLAS analyses do not extend to the region of the 
parameter space significant for us, that is m^ < 300 GeV and/or AM = m^—mLSP < m t- The 
reason is that signal efficiency relies in part on the presence of initial state radiation which 
is associated with large uncertainty. Therefore, our model is largely unconstrained by the 
LHC analyses, allowing for a stop as light as 100 GeV, the bound from LEP [41]. Dedicated 
searches for bi — > 6+LSP can potentially constrain the left-handed stop significantly. In this 
paper we will focus on the LHC sensitivity to the flavor mixed stop signals only. 

Beyond this particular topology, in general there is an insufficient coverage of the phe- 
nomenology of a light stop somewhat degenerate with the LSP in a large region of the available 
parameter space. Motivated by this lack of coverage several external analyses [42-46] have 
shown that LHC is sensitive to a broader region of the parameter space than the one explored 



by the public analyses. [42] and [43] analyze the LHC sensitivity for the stop four body decay, 
and a bound on its mass of around 250 GeV was estimated. Furthermore, [42] discuss the FV 
stop decay, focusing on the region m^ — m^ < my/ + mf, where the FV mode could dominate 
even in the MSSM. In [44] a model independent study was performed, and the FV stop decay 
into a light jet plus the LSP was investigated beyond this particular compressed region. A 
bound of ~ 300 GeV on the stop mass was derived from shape-based hadronic analyses and 
monojet searches [47, 48]. 

In this section we will study a simplified model where the stop decays with a 100% 
branching ratio to j'+LSP. We simulate the stop production in MadGraph5 [49]. The stops 
are then decayed and the decay products hadronized in Pythia [50], which also adds initial 
and final state showers. Since we are interested in hard jets arising from this radiation, we 
simulate t*t, t*t + j and t*t + 2j at the matrix element level and match it with the Pythia 
shower. We use PGS to simulate detector effects on our signal. The detector simulation and 
the QCD modeling are probably the strongest limiting factors of our simulation. We emulate 
relevant LHC analyses for our topology. 

4.1 Jets+MET searches 

Traditional jets+MET searches both of ATLAS and CMS are designed for QCD production 
for new physics particles such as squarks or gluinos, which subsequently decay into jets and 
the LSP - which escapes detection - leading to a missing energy signature. 

For spectra with a mass large splitting with the LSP, these searches are very effective. 
The limits on squark simplified models are crossing the TeV threshold. However, these limits 
are considerably weakened if only one quark eigenstate is light [51], even if the LSP is massless. 
The sensitivity of these searches is much worse for the compressed region (for instance, for 
our region of interest, rriq — ttilsp ^ mt)- Most of the data for these searches comes from 
high instantaneous luminosity samples. Consequently, the hadronic triggers relevant for these 
searches require a very high amount of hadronic activity (iJy + Et > 1200 GeV for relevant 
signal regions in the 4.7 ftr 1 ATLAS analysis, and H T > 500 GeV for the 4.98 ftr 1 CMS 
sample [52]). Therefore, these searches are very insensitive to our region of interest. The 
dedicated searches for the sbottom decaying into 6+LSP [53] could be potentially sensitive 
to our topology since a charm has a 10% probability to be mistagged as a 6-quark. It would 
be then be interesting to recast the dedicated searches for third generation squarks as done 
in [42] , but extending the region of the parameter space. The possibility of a large squark 
mixing with a c-quark final state could provide a motivation for an improved c-tagging. 

4.2 Shape-based searches 

Shape-based analyses like the CMS ay and the CMS razor analyses [54-56] are much more 
sensitive to our topology. Instead of only using the overall energy scale in the event, these 
searches also employ an event shape-based discriminant. This allows the use of shape-based 
triggers requiring much lower Ht thresholds. We derive limits in the stop-LSP mass plane for a 
large region of parameter space. We scan over the entire plane (beyond the compressed region) 



assuming a 100% branching ratio to jet + LSP. While this might be an oversimplification for 
a realistic model (for example, the single top channel, if open, might provide stronger limits), 
it allows us to investigate the efficiency of this particular topology in a model independent 
way. 

4.2.1 ax analysis 

The CMS SUSY hadronic searches employing ax [57] are sensitive to events with 2 or more 
energetic jets with missing energy. The ay variable was first discussed in [58]. Since the 
additional trigger on ax allows the Ht and px thresholds to be significantly lowered, these 
searches are more sensitive to our region over traditional jets+MET searches. 

In this type of analysis all events are clustered into a dijet topology containing two 
pseudo-jets, choosing the combination which minimizes the Ex difference between the two 
pseudo-jets. 

The variable ax is defined as, 

E% 



»t = ^ M T = V(^t + e t ) 2 ~ (p£ + A 2 ) 2 ~ (Py + Py) 2 (4- 1 ) 

where Ex for each pseudo-jet is defined as the scalar sum of px of its components. 

Jets considered in the analysis are required to have Ex > 50 GeV. Further, the two 
highest-Ii/r jets must satisfy Ex > 100 GeV each, and the highest Ex jet should be central, 
\i]\ < 2.5. Events with additional hard jets (Ex > 50 GeV) are vetoed if the jet \rj\ > 3.0. 
Events are also required to have significant hadronic activity Hx > 275 GeV. Finally, a 
requirement of ax > 0.55 brings the QCD backgrounds down to manageable levels. Events 
which pass these requirements are categorised according to their Hx- The px thresholds in 
the lowest two Hx bins are rescaled to keep the efficiency of the ax cut similar in each Hx 
region. 

The CMS analysis [57] does not cover the parameter space we are interested in. We 
emulate the ax analysis on our simulated sample (the limits are presented in figure 4 and 
figure 5). The experimental limit is calculated using a likelihood model to test for the presence 
of a variety of signal models. Here we will estimate the 95% upper limit on the cross section 
by the following formula, 

r, (AN-Lae) 2 , . 

X = — o = 3.84 (4.2) 

X N S M + 5 2 SM + Loe K } 

where e is the signal efficiency, L = 4.98 fb _1 is the luminosity, N$m is the expected number 
of events from SM and Ssm is the uncertainty in the expected number. The only relevant 
channel for this search is the b-quark jets channel. It turns out that a few bins have 
downward fluctuations in the background, and hence tend to lead to a tighter limit. To be 
conservative, we present both the observed and expected limits. We define AiV as, 

. , T expected limit , A „. 

AiV = ^ (4.3) 

Nobs — Nsm observed limit 
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Figure 4: Limits from the 7 TeV ax analysis in the (mr, m^sp) plane. Dashed red lines 
indicate observed limit, while the solid magenta line shows the expected limit. The stop is 
assumed to decay into jet and LSP with 100% branching ratio everywhere. 

Further, we do not attempt to combine limits from different bins, but consider a parameter 
point ruled out at 95% confidence if the production cross section for that point exceeds the 
limit from any individual bin. 

Our 95% level limits, computed for a single stop eigenstate decaying a 100% into jet and 
the LSP, are presented in the stop-LSP mass plane in figure 4. Our estimate excludes a single 
stop t below 250 GeV in the highly compressed region AM < mw + m& + niLSPi while for an 
increased mass splitting (AM ~ 150 GeV) the limit is around 350 GeV. Stop masses below 
~ 420 GeV are excluded for massless LSP. The interpretation of the actual limit within a 
realistic model with mixing is not straightforward since for AM > nit the decay mode into 
j'+LSP is no longer the dominant decay mode, and therefore the branching ratio suppression 
has to be taken into account. 

The left plot in figure 5 shows the 95% limit on the cross section for different mass 
splittings, 200 GeV < AM < 250 GeV. Also, we interpret these as limit on the branching 
ratio as a function of the stop mass. We notice that the cross section limits are relatively flat 
in mj, while the stop production cross section is steeply falling. Thus, the branching ratio 
suppression does not result in a significantly different limit. At the same time, when we add 
additional states at the same mass (for instance, in, ii, and hi all together or just ix and b£) 
the cross section gets multiplied by a factor of 2 or 3. By the same logic, these limits do not 
get much worse in this case. 

The plot in figure 6 shows the limits on the stop mass in the region AM > nit for a 
massless LSP for different values of the mixing angle. For a gravitino LSP, for mixing as low 
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Figure 5: Cross section (top) and branching ratio (bottom) limits from the 7 TeV (Xt analysis 
as a function of stop mass for various LSP masses, (solid black m^-— m^ 



10 GeV, dot-dashed 



blue: ?nr 
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150 GeV, dashed red: m^ 



GeV and m 3 / 2 = GeV). The gray band is 
the NLO stop production cross section at 7 TeV computed using Prospino 2.1 [59]. 

as 3% the entire stealth region mj < 250 GeV is excluded. For a massless bino instead, the 
stealth region is ruled out for a larger mixing angle sin 2 9 ~ 0.1. For a heavier stop, dedicated 
searches for t+LSP decay mode start to be sensitive setting a limit which depends on the 
amount of mixing. Stops of mass up to 500 GeV are excluded even in the maximal mixed 
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Figure 6: The limit on the branching ratio t —> jLSP as a function of the stop mass (black 
solid). The different curves represent branching ratios for t — > jG (dashed) and t — > jB (solid) 
considering sin 2 9 = 0.03 (red) and O.l(blue). The LSP is always assumed to be massless. 

scenario [33]. Therefore, the natural window for a massless LSP and a mixed stop could be 
almost entirely covered. 

4.2.2 Razor analysis 

The razor is an inclusive analysis which divides the signal and background into various boxes 
based on the final states. A global limit can then be derived from a likelihood analysis over 
every box. The power of the analysis lies in the fact that the backgrounds in each of the 
boxes are well-fit by a simple function, hence allowing for a relatively small error in the 
background extrapolation into the signal region. In our case, the signal only populates the 
HadBox (defined below), and hence all limits will be derived from this box only. 

The razor search is designed to look for pair production of new heavy states which yield 
high energy jets. All events are clustered into two megajets to reduce every final state into 
the dijet topology. The two razor variables employed, Mr and R are then calculated from 
these megajets. The Mr razor kinematic variable is given as, 



M R =[{\^\ + \^\) 2 -^+rt 



\213 



(4.4) 



and is invariant under longitudinal boosts. Mr approximates the energy scale in the heavy 
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particle decay and peaks at, 

2 2 

m~ — mi c P 
M A = —^ — (4.5) 

The i? variable is denned as 

R = M*/M R , (4.6) 

where MS < Ma is a transverse invariant mass, defined as, 

"1 



M^ 



K E^ ss (rt +Pt) ~ ET SS ■ (P>T +PT 



(4.7) 



QCD jets typically have a small value of R, and the QCD distribution peaks at R ~ 0, 
whereas the heavy particle production can produce higher values of R. 

In the compressed region, much of the jet activity arises from QCD radiation, and hence 
the above observations are not directly applicable. The Mr and R 2 distributions depend 
weakly on the stop mass, but are largely determined by the jet pt selection. Therefore, it is 
remarkable that the razor searches are effective in this region of parameter space. In fact, we 
find that the razor analysis limits go all the way to the degeneracy limit (m^ = rnLSp), where 
the stop decay products are not visible, and all hadronic activity arises purely from initial and 
final state radiation. It has been pointed out that in this region, the stop search is analogous 
to the dark matter search ([43, 60]). For the same reason, the monojet analysis designed for 
dark matter search sets significant limits on the stop masses for a very compressed spectrum 
(AM < 30 - 40 GeV) [42-45]. 

As mentioned above, the only box relevant for our simplified model is the HadBox. An 
event is put into the HadBox if it does not fall into any other Box (involving electron and/or 
muon final states), and satisfies M R > 400 GeV and 0.18 < R 2 < 0.5, [54]. 

The crucial feature of this search driving the limits in our case is the relatively low pr 
thresholds. All jets with pt > 40 GeV and |?7| < 3.0 are considered. There is a further 
requirement on the two highest-pT jets to have px > 60 GeV. 

The limits are derived according to the procedure outlined in [43] . From each box (which 
in our case is just one box), a posterior- probability density function is derived for the signal 
cross section, 

f°° f 1 (b + L<je) n e~ b ~ Lat 

P(a) = / db / de : lognormal(e|e, <5 e )lognormal(6|&, 6b) (4.8) 

Jo Jo n\ 

where b is the actual background yield, b is the expected background yield and 6b is the error 
in the background modeled using a lognormal distribution. The corresponding quantities for 
e denote the efficiencies. An error of 30% was chosen for the signal efficiency to account for 
detector simulation errors, a is the cross section, n the observed number of events and L is 
the available luminosity. The 95% limit is then straightforwardly derived from the cumulative 
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Figure 7: Limits on our scenario from the 7 TeV razor analysis [54] in the (mj, rriLSp) plane. 
The stop is assumed to decay into jet and LSP with 100% branching ratio everywhere. 

distribution of P(cr). In figure 7 and figure 8 we show the limits obtained. Our results are 
consistent with results in the literature [43, 44] to the extent they can be compared. 

We present the limits in the stop-LSP mass plane in figure 7, and the limits on the cross 
section and the branching ratio in figure 8. Similar conclusions as the ax search hold. The 
limit is insensitive to the stop mass, and hence branching ratio suppression or multiplicity 
enhancements do not change the limits significantly. 

Extending the razor 

The razor analysis appears to be applicable beyond its design topology. This was already 
noted in previous analyses [43, 60]. In particular, it is one of the strongest bounds on the 
signatures considered in this paper. We comment briefly on the possibility of extending the 
current analysis strategy to improve the sensitivity to our signal. 

In figure 9 we plot the two-dimensional histogram of the number of events expected at 
the LHC at 5 fb _1 as a function of the razor variables, rriR and B? for two mass spectra as 
examples. We see that for mj = 500 GeV, m ~ = GeV, the events populate high tur bins, 
and hence would tend to pass the analysis cuts. For m^ = 500 GeV, m~ = 300 GeV, the uir 
spectrum is significantly softer. However, the events are pushed towards larger E? values, 
where the QCD background is expected to be minimal. Therefore, extending the razor signal 
region to include the low m^ - high R 2 region could potentially improve the sensitivity to the 
signal. 
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[59]. 
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Figure 9: Two-dimensional histograms for the number of events at 5 fb 1 in rriR and R 2 
plane for m t - = 500 GeV, m^=0 GeV (top) and m t ~ = 500 GeV, m^ = 300 GeV (bottom). 

5 Conclusions 

We have explored a scenario where the third generation squarks mix significantly with either 
the first or the second generation. This arises generically (both in the left and the right- 
handed sector) in models with an approximate R-symmetry, but at the same time, these 
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results also apply to MSSM scenarios beyond MFV where the right handed stop mixes with 
the charm (for instance in the scenario described in [33]). These frameworks are especially 
interesting in light of the recent LHC null results for degenerate squarks. Also, independent 
of theoretical arguments, it is important to make sure that the LHC has a good coverage of 
SUSY scenarios beyond MSSM. 

In this paper we have focused in particular on a light stop somewhat degenerate with the 
LSP. We study the FV decay topology, t — > j LSP, which dominates in an extended region of 
the parameter space for moderately large mixing. This region is unconstrained in the public 
analyses of ATLAS and CMS. However, our emulation of the CMS razor and «t searches 
[54, 61] seem to be sensitive to this region of the parameter space and are able to rule out a 
significant portion of it. This is in agreement with recent work on light stop phenomenology, 
[42-45] . For a massive bino we find that the bounds on the stop mass varies between 250 GeV 
and 350 GeV depending on the mass splitting between the stop and the bino. For a massless 
bino the bound on FV decay mode rules out all the stealth region (rrif < 250 GeV) as long as 
mixing is around 10%. For a massless gravitino LSP, a stop below 250 GeV is ruled by the 
FV decay mode for an even smaller mixing, 3%. For heavier masses the dedicated searches 
for t—tt LSP start to be efficient, setting a bound in the 500-560 GeV range (depending on 
the mixing). 

These light stop bounds involve extending searches in regions beyond where they were 
designed to operate. There is a potential danger that this might not always be possible. For 
instance, the CMS razor analysis [62] attempting to put limits in the mass region rrig < 300 
GeV and vaq — m^sp < 150 GeV finds that this region has a sensitivity to the ISR modeling in 
simulation of signal events above a pre-defined tolerance; hence no interpretation is presented 
for these model points. This is an important piece of information which draws attention to the 
challenges for probing this region. We encourage the experimental collaborations to attempt 
to cover the light mass and/or somewhat compressed region for all the searches which are 
potentially sensitive to it such as for example the «t analysis and the razor analysis [54, 61]. 

It is important to extend the LHC searches for light stops to all possible channels 
unrestricted by theoretical bias. For instance, we show that a light stop can have signatures 
similar to a first or second generation squark in a significant region of parameter space. This 
is especially relevant for the 7-8 TeV run of the LHC, since at 13 TeV, the QCD backgrounds 
for jets+MET searches get even worse relative to the background for third generation final 
states. This could make our topology challenging to look for and potentially a blind spot for 
light stops. 
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